Multi-channel scattering and Feshbach resonances: 
Effective theory, phenomenology, and many-body effects 
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A low energy effective theory based on a microscopic multi-channel description of the atom-atom 
interaction is derived for the scattering of alkali atoms in different hyperfine states. This theory 
describes all scattering properties, including medium effects, in terms of the singlet and triplet scat- 
tering lengths and the range of the atom-atom potential and provides a link between a microscopic 
description of Feshbach scattering and more phenomenological approaches, ft permits the calcula- 
tion of medium effects on the resonance coming from the occupation of closed channel states. The 
examination of such effects are demonstrated to be of particular relevance to an experimentally 
important Feshbach resonance for 40 K atoms. We analyze a recent rethermalization rate experi- 
ment on 40 K and demonstrate that a measurement of the temperature dependence of this rate can 
determine the magnetic moment of the Feshbach molecule. Finally, the energy dependence of the 
Feshbach interaction is shown to introduce a negative effective range inversely proportional to the 
width of the resonance. Since our theory is based on a microscopic multi-channel picture, it allows 
the explicit calculation of corrections to commonly used approximations such as the neglect of the 
effective range and the treatment of the Feshbach molecule as a point boson. 
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I. INTRODUCTION 

The use of Feshbach resonances to vary the interac- 
tions between atoms obeying Fermi statistics has be- 
come a matter of interest in the field of ultracold trapped 
atoms 0. One reason is that these systems permit the 
experimental exploration of potentially strongly corre- 
lated Fermi systems in which the magnitude of the scat- 
tering length a is larger than the interparticle separation. 
The study of such strongly interacting Fermi systems is 
relevant to several fields including condensed matter and 
nuclear physics and astrophysics. 

It is therefore important to understand in some de- 
tail the physics of atoms interacting via a Feshbach res- 
onance. Theoretical studies of the two-body physics of 
Feshbach resonances are typically performed in the con- 
text of vacuum multi-channel scattering theory using a 
microscopic model of the atom-atom interaction |2( • On 
the other hand, studies of low energy many-body prop- 
erties often require the use of effective theories with a 
simplified form of the atom-atom interaction in order to 
make calculations tractable. A precise link between such 
effective theories and microscopic multi-channel models is 
therefore crucial in order to incorporate the correct two- 
body physics in many-body calculations. In the atomic 
gas community, effective theories are of two types. In one, 
the presence of closed channels is simply neglected, and 
only the resultant (resonant) scattering length is used as 
an input parameter. This approximation is adequate for 
the calculation of certain many-body properties close to 
broad resonances where the wave function of the Fesh- 
bach molecule is dominated by atom pairs in the open 
channel In the other, the presence of a two-body 



bound state in the closed channel (which leads to the 
resonance) is put in by hand in the form of a point bo- 
son 0,0- The parameters appearing in such theories are 
then the magnetic moment of the Feshbach molecule and 
the molecule-atom coupling. 

The treatment of the bound state in the closed chan- 
nels as a point boson is, however, incorrect in princi- 
ple since it is really a composite two-body object. The 
main purpose of the present paper is therefore to develop 
a systematic, improved low energy effective theory for 
the multi-channel scattering of alkali atoms in which the 
presence of a closed channel molecular Feshbach state is 
not put in by hand as a point boson but rather emerges 
naturally as a composite object consisting of two inter- 
acting fermionic atoms. We will show that such a theory 
permits tractable calculations of low energy many-body 
properties including medium effects due to the occupa- 
tion of states in the closed channels. The examination 
of such effects will be demonstrated to be of particu- 
lar importance for the experimentally well-studied reso- 
nance of 40 K atoms. It will be demonstrated that our 
approach can be reduced to the well-known Bose-Fermi 
model when the two-fermion nature of the bare molecule 
in the closed channel can be ignored. This theory can be 
formulated in two equivalent ways: One uses the micro- 
scopic parameters describing the atom-atom interaction 
as input; the other takes physical parameters describing 
the low energy scattering as input including the posi- 
tion and width of the resonance and the magnetic mo- 
ment of the Feshbach molecule (including the effects of 
dressing by the open channel). This provides an appeal- 
ing link between the microscopic parameters often used 
in two-body multi-channel calculations and the physical 
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parameters used in low energy effective many-body the- 
ories. We will use this effective theory to analyze recent 
rethermalization experiments on 40 K atoms and propose 
that further experimental study of rethermalization rates 
at higher temperatures can provide a determination of 
the magnetic moment of the Feshbach molecule. Finally, 
the energy dependent Feshbach interaction will be shown 
to introduce an effective range inversely proportional to 
the width of the resonance. This effective range can be 
neglected (and universal behavior close to resonance ex- 
pected) only in the case of broad resonances. 
Note that we set fi = 1 throughout the paper. 



II. LOW ENERGY EFFECTIVE THEORY 

In this section, we develop a low-energy effective theory 
for the scattering of alkali atoms in different hyperfine 
states. Consider a gas of alkali atoms in a given magnetic 
field B along the z-direction. Incoming and outgoing 
scattering channels are characterized by the eigenstates 
of the single particle hyperfine Hamiltonian 



H sp i n — A I 



S + CS,+DI. 



(1) 



where A is the hyperfine constant, C = 2{IbB, D = 
—/in B/I with n n the nuclear spin magnetic moment, and 
S and I are the electron and nuclear spin operators re- 
spectively. We denote the single particle eigenstates of 
the hyperfine Hamiltonian for a given B- field as \a) with 
pinl^) — € a |o;). The eigenstates can be labeled by their 
total spin, F = S + I, f = 7 ± |, and its z-projection, 



TTlF 



J + L e.g., 



\F, Trip). The main 



contribution to the atom-atom interaction driving the 
scattering process is the electrostatic central potential 



V(r) 



V.(r)+3V t (r) 



+ [V t (r)-V s (r)]S 1 -S 2 (2) 



where V s (r) and Vt(r) are the singlet and triplet poten- 
tials and Si and S2 are the spins of the valence elec- 
trons of the two alkali atoms |6J. We ignore here the 
small magnetostatic and magnetic dipole-dipole interac- 
tions between atoms. The potential J2J induces transi- 
tions among two-particle states with the same total z- 
projection of F . The properly anti-symmetrized two- 
particle state with the lowest energy constitutes the open 
channel \o) — \ai,at2). The states of higher energies 
which can be coupled to |o) by the potential @ form 

a set of closed channels |c^™') = lag™"*, a^). The thresh- 
old energies for the closed channels are then E { ^{B) = 



- e Ql , and they depend on the magnetic 



field B through H sp i n . 

As a practical illustration of our approach, we consider 
in the following the s-wave Feshbach resonance in a cold 
40 K gas at B ~ 201.6 G, which has frequently been 
studied experimentally, e.g., Refs. 0,0. The open chan- 
nel for this resonance has atoms in the states ||, — |) and 



||,— |)- The only other state in the ground-state man- 
ifold to which the open channel can couple through the 
dominant electrostatic potential has atoms in the states 



1 9 



||,-|) and ||,-|). Thus, \o) = % 
only closed channel is |c) = || — |, | 



|, I — |), and the 
- 1). Note that the 

1 O O \ . 

m com- 



2 > 



open and closed channels have the state 
mon. Thus, we can immediately conclude that any initial 
occupation of the state ||, — §) will have direct effects on 
the closed channel Feshbach state. Such effects cannot 
be described by the usual Bose-Fermi models which treat 
the Feshbach molecule as a point Boson. It was demon- 
strated in Ref. using a schematic mean field theory 
that these effects can be important for BCS-BEC cross- 
over physics. Here, we will demonstrate that the occupa- 
tion of the closed channel state can have significant effects 
on scattering properties using an effective theory based 
on a microscopic description of the atom-atom interac- 
tion. The threshold energy for this resonance is given 
by 



E: 



th 



Ei 



Ei__i_ =-\\A + V / 32A 2 + [2C -2D- 7Af 



Es._r = -i [A - y / 32A 2 + [2C -2D- 7Af 



(3) 
IAD], 
IAD]. 



With the values of A = -1.37 ■ 1(T 2 K, C = 1.34 
1CT 4 5/G K, and D = 1.19 • 1CT 8 5/G K from Ref. 
we can Taylor expand around Bq = 201.6 G to obtain 



E th ~ 0.084/f + 1.78/j, B (B - B a 



1 + 0.0188 



B-B 
100G 



(4) 



which is approximately linear in B — Bq for fields near 
the resonance. 

We emphasize that our effective theory is not limited to 
two-channel problems. It is readily generalized to more 
than two channels if appropriate for other atomic reso- 
nances. 

The total Hamiltonian describing a system of particles 
in the scattering states |o) and |c) interacting via V(r) is 



K,q 



Kc(q',q) Ko(q',q) 
V oc (q',q) K, (q',q) 



*q(K) (5) 



^K/2+qa 4 ^K/2-qa 3 ' ^K/2+qa 2 ^K/2-qai ] 



where *^(K) 

describes the creation of an atomic pair in the closed 
and open channels with center of mass momentum K 
and relative momentum q. Here, e a .k — ( a + k 2 /2m 
with m the mass of the atoms, and V is the volume of 
the system. The interaction matrix elements of V(r) in 
@ between the different scattering channels are denoted 
Vij (q', q) with i,j = o,c. They depend on B through the 
-B-dependence of the states \o) and |c). 

The scattering process is described by the scattering 
matrix, Tij(u>, K, q', q), which depends on the center-of- 
mass frequency and momentum (w, K) and the relative 
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momenta q and q' in the incoming and outgoing channels 
respectively. The scattering matrix obeys in the ladder 
approximation the Lippmann-Schwinger equation 



T T 

-L CC ± CO 

T T 

± OC -L OO 



v 00 



G c 
G n 



(6) 



where the propagator for two atoms in the open channel 
is 



G («,K,q) 
and 

G c (w,K,q,S) = 



l-/ Ql (±K-q)-/ Qa (±K- 



i X K 2 q 2 

uj + id — -r — 

4 m m 



l-/ a3 (|K-q)-/ a4 (|K + q) 
u + id-E th {B) - 



El 

4 m 



(8) 



is the propagator for a pair of atoms in the closed channel. 
Here f a is the Fermi distribution function appropriate for 
the hyperfme state \a). Note that © is a matrix equation 
both in the relative momenta (q', q) and in the channels 



at, and (implicitly) rc- Using QlOfl. one can readily calcu- 
late the low energy effective interaction matrix elements 
(i\U\j) for the scattering channels Note that these 
matrix elements depend on B through the independence 
of the eigenstates of QJ. For 40 K, a simple calculation 
expressing the hyperfme states |9/2, — 9/2), |9/2, — 7/2), 
and 1 7/2, —7/2) in terms of the eigenstates of the electron 
and nuclear spin for a given B-field yields 





'u cc 


Uco 


An 


(7) 




U 00 _ 


to 



c-?a s -{-at at — -a s 

Of — a B ccja s +at 



\/l+c 7 ^/l+c 



1 + C9 



(11) 



where c 7 = {2x - 7 - V81 - 28x + Ax 2 ) 2 / 32 and c 9 = 
(2x - 7 + V81 - 28x + 4x 2 ) 2 /32 with x = (C - D)/A 
incorporate the B dependence of the effective interaction. 

With (JTOJ [or (fTTfl for 40 K], we have obtained a low 
energy effective interaction. Since U is a solution to the 
scattering problem in vacuum, it includes the influence 
of the non-resonant coupling to high momentum states 
in both open and closed channels. 



B. Scattering matrix 



A. Effective interaction 

We now define an effective interaction U appropriate 
for calculating the low energy scattering properties of 
the system. This is done by considering the Lippmann- 
Schwinger equation in vacuum for lo = and ignoring 
the hyperfme splitting of the atomic levels 



(9) 



~u cc 


U co ' 


-1 


'Vac 


Vco 


-1 


■ G vac q 




Uoo_ 




Voc 


Voo_ 




G vac 



where G vac ((?) = (iS — q 2 /to) -1 is the vacuum pair prop- 
agator ignoring the hyperfme splitting. When the scat- 
tering is dominated by the presence of a single Feshbach 
state in a closed channel, the effective interaction medi- 
ated by the closed channel state is rank one separable: 
Its Fourier transform can be expressed in terms of the 
product of the Fourier transforms (4>k\V\4> m ) , where |</> m ) 
is the molecular state and \4>k) is a plane wave state of 
the two atoms 0- For the low energies relevant for di- 
lute atomic gases, the typical wavelength q^ 1 of the open 
channel scattering atoms is such that qrc <SC 1 where rc 
is the characteristic length scale of (r|V"|(/> m ). For these 
long wavelengths, the solution of @ can be written from 
© as 



#(q',q) 



4tt 



3ot 



+ (a t - <z s )Si • S 2 



9(q')9(q) 
. ( 10 ) 

where a s and at are the scattering lengths for the singlet 
V s (r) and triplet V t (r) potentials, respectively. Cut-off 
effects of the potential (for rc > 0) are contained in the 
form factor g(q) — > for qrc — * oo. The atom-atom po- 
tential is thus characterized by the three parameters: a s , 



T T 

± CC -L CO 


-1 


'Ucc 


Uco 


-1 


AG C 





T T 

± OC -L OO 




Uoc 


Uoo_ 







AG 



We can now use U as an effective low energy interaction 
for the exact solution of the scattering problem. The 
resonant scattering matrix can be expressed in terms of 
U using © and © with the result that 



(12) 



where AG„(w,K,q) = G Q (w,K,q) - G vac (q) and 
AG c (w,K,q,B) = G c (w,K,q,B) - G vac (<j). The pair 
propagators AG G and AG C provide corrections to the 
scattering matrix arising from the hyperfme level struc- 
ture which leads to the presence of a Feshbach resonance 
and corrections due to the presence of a medium in both 
channels. Note that since the effective interaction U in- 
cludes the coupling to high energy states, our theory as 
defined in l|12|l is properly renormalized in the sense that 
it yields finite results even for rc — corresponding to a 
5-function interaction. 

With U given by (| 1 0|) . we can readily solve l|12|l to 
obtain the scattering matrix in the open channel 



Of 

where 

n (w,K) = 

and 



U 



i - u cc n c 



(13) 



d 3 q 

(27T) 



[G (u;,K, q )-G vac (q)}g 2 (q) (14) 



n c (a;,K, J B) = J ^ [G c ( w ,K,q,B) - Gr™(a)}g 2 (q). 

(15) 
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The integrals in ^14 (1 - l|15 fl can performed straightfor- 
wardly in the vacuum case for a zero range (rc = 
0) potential with the result that n™ c (w,K) = 
— im 3 / 2 v / w/ (47r) + 0{\fmuj rc) where Co = uo — K 2 /(4m) 
is the energy in the center of mass frame. Similarly, 
W c ac {uo,K,B) = n™> - E th (B),K] +0{V^E^r c ). 



C. Feshbach resonance 

Consider now the case when a Feshbach resonance is 
present at B = Bq. For fields close to Bo, it is standard to 
parameterize the zero energy vacuum scattering matrix 
as 



rpvac 
oo 



AB 

B-B 



(16) 



where a^ g is the (non-resonant) background scattering 
length and AB the width of the resonance. We see from 
(113(1 that the zero energy (lo = 0) scattering matrix in a 
vacuum is 



T vac (0) = [/ oo + |[/ oc | 2 



n™ c (o) 



i-c/ cc n™ c (o) 



(17) 



since n™ c (w = 0) = 0. A comparison of (T^ with itTTI) 
allows us to express the phenomenological low energy 
scattering parameters etb g , AB, and Bq in terms of the 
microscopic parameters a s , a t and rc which character- 
ize the atom-atom interaction. The resonant part of the 
scattering coming from the coupling between the open 
and closed channels is clearly given by the second term 
in (|17|l . which depends on the B field through U and 
n™ c (0) . The position of the resonance Bq is thus deter- 
mined by 

l-[/ cc (S )nr c (5 ,w = 0) = 0. (18) 
A straightforward pole expansion in B — Bq of l|17|l yields 



The = 



47rab p 



m 



2U 0C d B U c - U: 



U 2 

II - — 

i d 2 B (u c<l -u 2 nm 



d B (u cc -u? c n™) 

for the background scattering and 



(19) 



j=0,B=B o 



T hs AB = - 



Uqc(Bo) 2 /U cc (Bq) 
d B [l-U cc {B)^{BM\i 



(20) 



for the width of the resonance. With {Hty. (fTT^fl - ftfy. 
(119(1 . and 1(20(1 . we have finally arrived at our low energy 
(kre <C 1) effective theory. Equation (|T3|) gives the full 
low energy/momentum dependent multi-channel scatter- 
ing matrix as a function of the parameters a s and a t 
(and rc) through U in l(10|) . Medium effects due to the 
occupation states in both open and closed channels are 



included explicitly included through n o and n c . Further- 
more, 1(18(1 . ((19(1 . and ((20(1 express the physical parame- 
ters dbg, Bq, and AB describing zero energy scattering 
in vacuum in terms of a s and at and rc- Alternatively, 
these equations can be used to fix a s and at and rc from 
experimental data. 

We emphasize that our theory is finite and well defined 
even for rc — (i.e., g(q) = 1). The form factors g(q) 
have been retained in order to describe finite range ef- 
fects. These can important since the hyperfine splitting 
energy E t \ Y can be comparable to h /(mr c ). 



III. MAGNETIC MOMENT OF THE 
FESHBACH MOLECULE 

In the spirit of Landau Fermi theory, it is often 
convenient to express the full energy- and momentum- 
dependent multi-channel scattering matrix in terms of 
physical observables only. To do this, it is necessary to 
introduce one more parameter: the magnetic moment of 
the Feshbach molecule in a vacuum. The energy of the 
molecule in a vacuum is determined by the poles of the 
full pair propagator in the closed channel: 



d£(k,lj) = 



IT 



II TT vac 



D 



bare 



u 



IT 



1 



TJ TTvac 

U OoTT 



(21) 



The first term, = ILY ac_1 — U cc , describes the prop- 

agation of a pair of atoms interacting in the closed chan- 
nel only (i.e., the "bare" molecule) while the final term 
describes the coupling to the open channel ("dressing"). 
The molecule propagator is illustrated diagrammatically 
in Fig. ^ To obtain an expression for the magnetic mo- 
ment of the molecule for fields close to Bq, we expand 
((21(1 about B — Bq and Co — to obtain 



Afx(B - Bq) 

ABA^l 



(22) 



where Co = lo — K 2 /Am is the energy in the center of mass 
frame as before and 

-l 



Afj, 



u cc }\ 



-i 



= zA^ba 



(23) 



is the magnetic moment of the molecule. We have defined 



An 



bare 



d B [w c 



cc\ I B 



d B E th + 2E th U- 1 d B U c , 



(24) 



using W c ac {B, to) cx [E th (B) - lo} 1 / 2 for r c = to calcu- 
late the partial derivatives. In addition, 

z- 1 = i-[/ oc 2 a u AnrRnr^ 

= 1 - TbgABA^barcCLAIir , (25) 
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(|14|l to calculate Alio a s 



Dressed Molecule: 
&&&&&& = -/\/wv + 



Bare Molecule in closed channel: 




\U + 

■ cc 

-i» 



FIG. 1: The full (dressed) molecule given by 12H is indi- 
cated by the double wavy lines whereas the bare molecule in 
the closed channel, excluding coupling to the open channel, 
is given by the single wavy line. Fermions are indicated by 
straight lines, and the intrachannel couplings U 00 and U cc are 
indicated by dashed lines. The coupling U oc between the open 
and closed channels is indicated by •. 



where we ha ve use d T hg AB A^w^n™ " 1 = U 2 oc fol- 
lowing from (|20I24|I and explicitly separated the thresh- 
old \JTo term in n o by writing 



7 3 / 2 



-in 



An (w). 



(26) 



The term AII contains the corrections to the pair prop- 
agator in the open channel coming from the finite range 
(rc 7^ 0) of the potential. Equation H23fl - I|25|> provides a 
clear physical interpretation of the various contributions 
to the magnetic moment of the molecule. Equation 124|) 
gives the magnetic moment of the bare molecule when 
the interaction with atom pairs in the open channel is 
ignored. It consists of a contribution from the magnetic 
moment difference between the open and closed channels, 
dsEth and a contribution from the magnetic field depen- 
dence of the interaction matrix elements in the closed 
channel. For the 40 K resonance discussed above, we have 
from (@J d B E t h = 1.78/j b . The second term in can 
be estimated as 



2E th U- c 1 d B U cc ^0A2^ B 



a i 



0.89- 



Oi 



The factor z given by l|25(l gives the renormalization of 
the magnetic moment due to coupling of the molecule to 
high momentum atom pairs in the open channel (exclud- 
ing the threshold effect given explicitly by the -^/tj-term) . 
This renormalization was first included systematically in 
a low energy effective theory in Ref. |5j where it was ar- 
gued that it can, in general, be significant. To see this, 
we use a simple cut-off form factor g(q) = 0(r ( ^ 1 — q) in 



2 . 

-A/i b 

7T 



1-0.35 



A/i b 



ABmab g rc 
are AB m a h „ r c 



10G 40u 100A 10A 



(27) 



where 1 u is the atomic mass unit. For the Feshbach res- 
onance of 40 K considered in this paper, AB ~ O(10G) 
and a bg ~ O(100A) (see next section). It then follows 
from Ij27(l that the renormalization of the magnetic mo- 
ment of the molecule due to the coupling to high momen- 
tum atom pairs in the open channel can be significant 
(i.e. z < 1) if rc ^ 10A. Likewise, the renormalization 
can be significant for the very broad resonance for 6 Li at 
B ~ 830G with AB ~ 300G and a bg ~ O(1000A) \V2}. 
However, for the narrow 6 Li resonance at Bo ~ 543G 
and AB = 0.23 G studied in Ref. 11], |23 indicates 
that the renormalization of the magnetic moment away 
from threshold described by z can be ignored, i.e. z ~ 1. 

Having introduced the magnetic moment of the 
molecule, we can now express the full frequency- and 
momentum-dependent scattering matrix in a medium in 
terms of the physical parameters ab g , Bq, AB, and A/i. 
To do this, we rewrite the open-open channel scattering 
matrix given by (|13fl in the form 



T — 

± on 



Uo 



Ul/[1 - U 00 IL } 2 



u cc - u, 



n„ 



i - u 00 n 

Uoo 

1 - Uoollo 1 - UooUo ~ 1 - Uoollo 



-D- 



oc i-u 00 n 

Uco 



(28) 



The first term describes the scattering via the open chan- 
nel interaction only, and the second term describes scat- 
tering via the closed channel. The in-medium molecule 
propagator D is given by l|21|) with the replacement 
rjvac _^ jj This splitting of the scattering into open 
channel and closed channel processes is illustrated dia- 
grammatically in Fig. [5] Comparison with l|21|l reveals 
that the second term in l|28ll describes the presence of the 
Feshbach molecule (including medium effects) leading to 
the resonance. Note that both terms in (|28|l contribute 
to the background scattering given by (|19fl . Indeed, a 
Taylor expansion around u> = and B — Bo and com- 
parison with the vacuum expression for the scattering 
matrix (|16f) yields 



T — 

± no 



Q+h(ui)-Aij.(B-B ) 



(29) 



- r bg n«rf H 



where Hg(u), K) is the pair propagator in the open chan- 
nel assuming rc — 0. The term 



U c (u,K,B) -n^ c (ff,0) - d u II?*(B ,u)\ u> 



n c (w,K,s) -n™ c (s,o) 

d u U^(B ,u)\o 



(30) 
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\u + 


\u 00 \ 


i oo 






FIG. 2: The open-open channel scattering matrix split into 
scattering in the open channel only (first line) and scattering 
via the closed channel (second line) as given by 1281 . 



in (I29fl describes effects coming from the fact that the 
Feshbach state is a composite two-fermion object includ- 
ing medium effects in the closed channel. 

With (|29[1 . we have expressed the full scattering ma- 
trix in terms of the physical parameters etb g , Bq, AB, 
and A/i. These, in turn, can be calculated from a s and 
a t with the aid of JTBJ-lGnjl and Thus, (HU and 

(|13[l are two formulations of the same result, the for- 
mer expresses the in-medium scattering matrix in terms 
of the microscopic parameters characterizing the atom- 
atom interaction; the latter provides a formulation valid 
close to the resonance in terms of the physical observ- 
ables characterizing the resonance. Our effective theory 
thus provides a link between these descriptions. 



IV. MEDIUM EFFECTS 

We now consider the shift of the resonance position 
due to the presence of a medium. To do this, we rewrite 
(EH as 



+ - 



1 - Tbglljff 
T hg AfiAB/[l-T hg Il^] 2 



Co + h{uo) - A/j,(B - B a ) - TbgAfiAB- 



(31) 



T bE n" 



We see clearly that the term T hg AnABU^ { in the de- 
nominator describes the effects of the presence of atoms 
in the open channel whereas h(ui) describes the effects of 
the presence of atoms in the closed channel. 

We first calculate medium effects due to the presence 
of open channel atoms. We see from (|14|l that medium 
effects in the open channel shift 11 by 



d 3 q /ai(^K-q) + / a2 (^K + q) _ mk F 



(32) 



for K = and u> = where kp is the Fermi momentum 
for both open channel hyperfine states (we have assumed 
equal population) and ksT <C k F /2m. The resulting 
shift SB in the position of the center of mass K = 
resonance from its vacuum position Bq is then readily 
determined from (|31|l 



SB = k 

7T 



F^bg 



AB. 



(33) 



Taking a density of n = 2 x 10 12 cm~ 3 , a^g = 290a , 
and AB — 8G relevant for the 40 K resonance considered 
here, l|33|l yields SB — — 0.61G for the shift of the posi- 
tion of the resonance due to the occupation of states in 
the open channel. Of course, in a typical radio-frequency 
(RF) experiment which probes the scattering of pairs of 
atoms within the Fermi sea with varying center of mass 
momenta K , it is necessary to average (1141) and l|31() 
appropriately over momenta in order to obtain the aver- 
age medium shift of the resonance position. 

The effects due to the population of the closed channel 
states can likewise be calculated from h{ui). The medium 
shift of h(u>) is obtained from l|3U|l and (|15|l as 



d 3 q / a3 (|K- q) + / a4 (|K + q) 
(2?r) 3 E th 



K' 2 
4 m 



n c 

Eth ' 



(34) 



where n c is the density of closed channel atoms and 
E th > (37r 2 n c ) 2 / 3 /2m by assumption. The shift of the 
resonance position due to the occupation of closed chan- 
nel states follows from l|31ll as 



SB 



8tt 



2n c U c , 

A/ibar. 



(35) 



For the 40 K resonance, an assumed total density of n = 
2 x 10 12 cm~ 3 corresponds to n c = 10 12 cm -3 since the 
open and the closed channels share one hyperfine state. 
Using E t h = 0.084if from and assuming A/ibarc ~ Ms 
yields SB = — 1.7 x 1CU 6 G for the shift of the position of 
the resonance due to the occupation of states in the closed 
channel. The shift SB is very small because E t h ^> e_F 
for this resonance where tp is the Fermi energy. 

The medium effects are treated within the ladder ap- 
proximation in the present paper. Close to resonance 
when interaction effects are strong, it is not obvious 
that this approximation is adequate. However, the 
good agreement between the single channel thermody- 
namic calculations based on the ladder approximation 
and Monte Carlo results indicate that the ladder approx- 
imation includes the most important physics, even close 
to resonance [T^. 



V. COMPARISON WITH BOSE-FERMI 
EFFECTIVE THEORIES 

A comparison with other recent low energy effective 
theories for atomic Feshbach resonances is now in order. 
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Most of these theories are based on a Bose-Fermi model 
where the presence of the Feshbach molecule is put in by 
hand as a point boson 0, 0| . In the present approach, 
on the other hand, the Feshbach molecule emerges dy- 
namically and is treated correctly as a composite two- 
fermion object consisting of atom pairs in both the open 
and the closed channels. This allows for the descrip- 
tion of effects coming from the two-fermion nature of the 
molecule through II C in (fH?|) or the term h(oj) in 
(|30[l . The simpler Bose-Fermi theory, which describes the 
Feshbach state as a point boson whose propagator has 
linear (^-dependence, is recovered only when h(u>) = 0. 
Furthermore, the present theory provides expressions for 
all low energy scattering and bound state properties, in- 
cluding medium effects in both open and closed channels, 
through the parameters et s , at, and tq characterizing the 
atom-atom interaction and thus provides a nice link be- 
tween a microscopic (coupled-channels) description of the 
scattering and more phcnomcnological descriptions. This 
includes a physically clear picture of the various contribu- 
tions to the magnetic moment of the Feshbach molecule. 
Once a Sl a t , and tq are fixed, our approach affords a 
systematic description of Feshbach resonances in all the 
relevant channels. 



VI. RETHERMALIZATION RATE 

Finally, we analyze a recent experiment on the Fesh- 
bach resonance at Bq ~ 201.6 G for 40 K atoms in the hy- 
perfine states |9/2, —9/2) and |9/2, —7/2). In refs. 0,||, 
the elastic collision rate for 40 K atoms was measured near 
this Feshbach resonance using the following rethermal- 
ization technique: The gas was "heated" preferentially 
in one spatial direction, and the relaxation toward equi- 
librium with a uniform temperature in all directions was 
then followed in the time evolution of the rms cloud radii. 
The rethermalization rate was then extracted from an ex- 
ponential fit of the aspect ratio as a function of time. 

A variational expression for the relaxation time of tem- 
perature anisotropies can be given [|J as 



1 



(($ T r[$ T ])) 



(36) 



Here. 



«•••» = 



d 3 r 



d 3 p 



./o(l-/o) (37) 



where /o is the equilibrium distribution function for 
|9/2, —9/2) and |9/2, —7/2) atoms and where we treat 
the trap potential T4 rap (r) using the Thomas-Fermi ap- 
proximation with a spatially dependent chemical poten- 
tial /i(r) = [i- Vtrap(r). The populations of |9/2, —9/2) 
and 1 9/2,— 7/2) atoms are assumed to be equal. The 
deviation function appropriate for thermal anisotropy is 



4>7 



2^2 



■Y z -p i /3+- s -fr4z 



>i) 



(38) 



where lo z and u>± are the trapping frequencies in the axial 
and transverse directions (we assume an axially symmet- 
ric trap). The linearized collision operator is 



+ S'J/O/ftl - /°')(1 - /°'), 



Vl 



(39) 



where da/dQ, is the on-shell (ui — p 2 /2m + pi /2m) differ- 
ential cross section and fl is the solid angle for the direc- 
tion of the relative outgoing momentum p'. = (p' — p'i)/2 
with respect to the relative incoming momentum p r = 
(p — Pi)/2 6j. The calculation of the integrals in (I36*|) 
proceeds as in Ref. [l^. In the classical limit, these inte- 
grals can be evaluated analytically yielding 



1 

tt 



4 

5N 



1/2 



d 3 rn(r) 5 



(40) 



where N is the total number of trapped atoms (in both 
hyperfme states), n the total density of atoms, and 



dxx 7 a(\/ mksTx)e x 



(41) 



The factor of x 7 in l|41|) comes from the weight function 
$t appropriate for temperature relaxation. Note that 
(|40|l differs from the averaged cross section arising from 
the collision relaxation time. The latter results from a 
thermal average of the cross section weighted with the 
relative velocity only, which corresponds to the factor x 3 
in gTJ. 

We now use a = m 2 \T 00 \ 2 /4ir with T ao given by lf2*9"|) 
to analyze the experimental data. In Refs. 0,0, the ex- 
periments are performed at temperatures T ~ 2Tp and 
densities n ~ 10 12 cm -3 . We can therefore use the clas- 
sical results (|40f) and l|41|l . Likewise, T OQ in 129fl can be 
approximated by the vacuum on-shell scattering matrix 



T — 

-L on — 



4 IT „ 



1 



Afj.AB 



q 2 /m-&n{B-B ] 



(42) 



ia hg q 



where q is the relative momentum of the two scatter- 
ing atoms. At present, knowledge of the experimental 
value of the magnetic moment of the Feshbach molecule 
is limited. We therefore consider two cases of a large and 
a small value of A/i. We take A/j, = dsEth — 1-78 /is 
from Q as a large value and A/i = A/io = 0.118 /ib from 
Ref. |5j as a small value. Using l|41|l - H42|) . we can now an- 
alyze the experimental results |7J,|8j . For the two values of 
A/x above, we tune the parameters AB and a^ g keeping 
the resonance position fixed at B = 201.6 G. The results 
of our fits are presented in Fig. El For A/i = 0.118 \ib we 
find 

AB = 6.5 G and a hg = 286.7 a a . (43) 
The fit with A/i = 1.78 /j,b gives 

AB = 8.1 G and a bg = 296.7 a . (44) 
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FIG. 3: Fit to the experimental rethermalization cross sec- 
tion H with I4H - 1421 for two values of the molecule magnetic 
moment at T — 4.4 fiK. The fit parameters are given in l|4M[l 
for Afi = 0.118 hb and JUJ for Afi = 1.78 fi B - 
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FIG. 4: The rethermalization cross section l|41|l with l|42[l for 
two values of the molecule magnetic moment and the temper- 
atures T — 8 fiK and 20 fiK. The fit parameters are given in 
for Aft = 0.118 fi B and (£Q for Afi = 1.78 (la- 



in Ref. 0, the authors give the fitted values of the 
singlet and triplet scattering lengths as a s w 104.8 ao 
and at ~ 174 ao- Substituting these values in i|19fl. we 
find abg — 170 an, which is close to the value 164 ao re- 
ported in Ref. (l5|. This value differs somewhat from 
ours (abg — 290a ). The reason for this discrepancy can 
be the experimental uncertainty in determining the num- 
ber of atoms trapped, and/or the fact that the authors in 
Ref. Q use a different thermal averaging procedure than 
that given bv l|IT Jl [T ^ . 

The fit parameters (I43|) - (I44H do not depend strongly 
on the value of the molecule magnetic moment at the 
temperature T — 4.4 /iK. This is because AfiAB ^> k B T 
for both values of Afi: AfiAB j k B = 55fiK for Afi = 
0.118 fi B and A/j,AB/k B = 980fj.K for Apt = 1.78 fj, B . In 
fact, the overall shape of the rethermalization curve as 
a function of B — Bq (e.g. the position of the resonance 
and the magnitude of the maximum) becomes sensitive 
to A/j only when 

AfiAB - k B T. (45) 

To illustrate this, we show in Fig.0]the cross section l|41() 
evaluated for temperatures T = 8 [iK and T — 20 /iK and 
for the two values of A/Lt above. We see a clear difference 
in the shape of the curve depending on the temperature. 

The magnitude of the minimum of the rethermal- 
ization rate located at B — Bq ~ AS (for k B T <C 
A/iAB) is however sensitive to the temperature even for 
k B T A^iAB. This can easily be understood as fol- 
lows. If one ignores the q 2 /m term in (|42|l . the cross 
section is identically zero for B — Bq — AB for all q 
leading to a vanishing rethermalization rate. But from 




An [nj 

FIG. 5: The position B — Bo — Bmi n and magnitude (ffmln) 
of the minimal value of the rethermalization cross section as 
a function of the molecule magnetic moment Afi for T = 
4.4 fjK, T = 8 /iK, and T = 20 fjK. 

q 2 /m ~ k B T is follows from l|42|l that the minimum of 
the rethermalization rate at B — Bq ~ AB scales as 
k B T/ AfiAB for k B T < AfiAB. Thus, the minimum 
of the rethermalization rate is sensitive to the tempera- 
ture or the value of the molecule magnetic moment even 
for k B T < AfiAB . To illustrate this, we plot in Fig. 
the position £> m i n = B — Bq and the magnitude (er m i n ) of 
the minimal value of the rethermalization cross section 
as a function of the molecule magnetic moment Afi for 
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T = 4.4 fiK, T = 8 /iK, and T = 20 /zK. We see that the 
position of the minimum does not depend on A[i when 
k B T <C A/iAB whereas the magnitude of the minimum 
value of the rethcrmalization rate does. 

The temperature dependence of the rethermalization 
rate could be used to determine the magnetic moment 
of the Feshbach molecule: One can measure this mag- 
netic moment by fitting the measured rethermalization 
cross section as a function of temperature using the 
theory described above. A scan of temperatures from 
k B T < A^iAB to k B T > AfiAB is likely to yield a fairly 
accurate value of A/i. For k B T <C A^xAB, the Feshbach 
resonance gives rise to a sharp maximum and minimum 
in the rethermalization scattering cross section whereas 
for k B T 3> A/iAB, the rethermalization cross section is 
a much smoother function of the B field. Furthermore, 
if one is able to measure the magnitude of the minimum 
rethcrmalization rate which is sensitive to the tempera- 
ture for all T, the size of the magnetic moment can even 
be determined for k B T <C A/iAB. 



VII. EFFECTIVE RANGE 

We now demonstrate that the energy-dependent, 
Feshbach-mediated interaction introduces a new length 
scale: the effective range r e g. The effective range expan- 
sion for the vacuum scattering matrix is given by [Tt| 



T — 

-L on — 



4:770, 



m 1 



\q 2 ar eS 



i qa 



(46) 



where a = ab g [l— ABj (B— Bo)] is the resonant scattering 
length. This expansion can be obtained from l|42(l in the 
limit \B ~ B \ <C AB with the effective range given as 



mA/iA_Babg 



(47) 



The q 2 /m term in (14211 coming from the momentum de- 
pendence of the Feshbach molecule energy can therefore 
introduces a new length scale: An effective range given 
by H47fl . Equation l|47|l was also derived in Ref. [18( in 
the context of a Bose-Fermi model. 

Note that the effective range given by l|T7|) is nega- 
tive as a direct consequence of the energy dependence of 
the multi-channel interaction. In a recent paper, finite 
range effects on the physics of cold atomic gases in the 
BEC-BCS crossover regime was examined using a single 
channel model |l9j . The atom-atom potential used in 
Ref. |l9j has a positive effective range and thus does not 
provide a correct description of finite range effects for 
resonant multi-channel atom-atom scattering. Similarly, 
the effective range of the single channel interaction used 
in a recent Monte-Carlo study also has r c s > [2J}. In 
general, it is often difficult and sometimes impossible to 
obtain a negative effective range in single channel models 
of resonant scattering. This feature emerges naturally in 
multi-channel models. 



It is common to neglect the effective range term by 
adopting a simple approximate form of the unitarized 
scattering matrix 



T 



4na 



1 



m 1 + iqa 



(48) 



This approximation leads to the idea of universality. 
When \a\ — > oo, the only remaining length scale in the 
problem is the interparticle spacing, and the properties 
of the gas should therefore be independent of the micro- 
scopic details of the atom-atom interaction. The physics 
of interacting gases in the universal limit has been exam- 
ined by a number of authors [20 l I2 I I2ll2 ^| . However, 
we have now demonstrated explicitly that there is an ad- 
ditional length scale in the problem given by l|47|l . This 
length scale can be neglected only when calculating prop- 
erties with a typical wave length q for which qr c ^ -C 1; 
the thermodynamic properties of the gas are independent 
of r e ff only when 



(49) 



where r*i2 is the average interparticle spacing. For low 
temperatures where q ~ kp, using (|47|) . this condition 
can be written as 



I A/iAB I > 



fcf|abg| 



(50) 



The simple unitarized form for the scattering given by 
H48J1 therefore holds only for broad resonances as defined 
by (|5Ufl . In particular, one can expect universal behavior 
only for such broad resonances [13 ■ 

For the Feshbach resonance for 40 K at B ~ 201. 6G, 
we find AB ~ 10G and ab g ~ 290ao as explained 
above. Assuming Afj, ~ fi B and ep = 10~ 6 K, this yields 
A/iAB 3> ef/fep-Obg and one would expect the gas to ex- 
hibit universal behavior close to resonance. This is also 
true for the very broad resonance for 6 Li at B ~ 830G 
with AB ~ 300G. However, for the narrow 6 Li resonance 
at B a ~ 543G, AB = 0.23G, and a bg - 80a [13, we find 
tp ~ 3kpat,gfiBAB for ep = 10 _6 K. One can therefore 
not expect the effective range to be an irrelevant length 
scale for this narrow resonance. 



VIII. CONCLUSIONS 

A low energy effective theory for the in-medium scat- 
tering of alkali atoms was obtained using a microscopic 
multi-channel description of the atom-atom interaction. 
Expressions for all scattering properties in terms of a few 
parameters characterizing the atom-atom potential were 
derived thereby providing a link between a microscopic 
multi-channel description of Feshbach resonances and 
more phenomenological approaches. The presence 
of a Feshbach resonance emerged naturally from the 
theory as a two-particle bound state, and the energy 
dependence of the Feshbach-mediated interaction was 
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shown to introduce a negative effective range inversely 
proportional to the width of the resonance. Furthermore, 
medium effects arising from the occupation of both open 
and closed channel states were calculated. The theory 
was shown to allow the explicit calculation of correc- 
tions to commonly used approximations such as the 
neglect of the effective range and the treatment of the 
Feshbach molecule as a point boson. We also analyzed 
a recent rethermalization experiment on 40 K atoms and 
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